Localization and fluctuations of local spectral density on tree-like structures with 
large connectivity: Application to the quasiparticle line shape in quantum dots 
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We study fluctuations of the local density of states (LDOS) 
on a tree-like lattice with large branching number m. The av- 
erage form of the local spectral function (at given value of the 
random potential in the observation point) shows a crossover 
from the Lorentzian to semicircular form at a ~ 1/m, where 
a = (V/W) 2 , V is the typical value of the hopping matrix 
element, and W is the width of the distribution of random 
site energies. For a > 1/m 2 the LDOS fluctuations (with re- 
spect to this average form) are weak. In the opposite case, 
a < 1/m 2 , the fluctuations get strong and the average LDOS 
ceases to be representative, which is related to the existence 
of the Anderson transition at a c ~ l/(m 2 log 2 m). On the 
localized side of the transition the spectrum is discrete, and 
LDOS is given by a set of J-like peaks. The effective num- 
ber of components in this regime is given by 1/P, with P 
being the inverse participation ratio. It is shown that P 
has in the transition point a limiting value P c close to unity, 
1 — P c ~ 1/logm, so that the system undergoes a transition 
directly from the deeply localized to extended phase. On the 
side of delocalized states, the peaks in LDOS get broadened, 
with a width ~ exp{— const log m[(a — a c )/a c ] -1 ' 2 } being 
exponentially small near the transition point. We discuss ap- 
plication of our results to the problem of the quasiparticle 
line shape in a finite Fermi system, as suggested recently by 
Altshuler, Gefen, Kamenev, and Levitov. 



I. INTRODUCTION. 

The recent experimental observation of single-level 
peaks in excitation spectra of quantum dots plpj has mo- 
tivated theoretical interest in the problem of the quasi- 
particle life time induced by the Coulomb interaction 
in mesoscopic samples [J2|j3j]. It was found that the sin- 
gle particle levels with excitation energy E < Et have 
widths less than the one-particle mean level spacing A, 
and thus can be resolved. Here Et is the Thouless en- 
ergy, which has the physical meaning of the inverse time 
of spreading of a wave packet over the system, the ratio 
g = Et/A being the dimensionless conductance of the 
dot. Thus, the theory predicts that approximately the 
first g single-particle levels can be resolved, in agreement 
with experimental findings |l|,|2|,|4| . In a more recent pa- 
per B, Altshuler, Gefen, Kamenev and Levitov (AGKL) 
mapped the problem onto a tree-like tight-binding model 
in the Fock space. This model is very close to the Ander- 
son tight-binding model on the Bethe lattice, which was 
found [(Hi to undergo the localization transition. This 
allowed AGKL to conclude that there is a localization 



threshold E c (A < E c < Et) in the problem, so that 
the states well below E c are strongly localized, i.e. given 
by the bare single-particle states with small perturbative 
admixture of many-particle states. 

The qualitative arguments used by AGKL are not suf- 
ficient, however, for a more precise determination of the 
behavior of the spectral shape of the quasiparticle peak in 
the vicinity of the localization threshold. On the other 
hand, the critical behavior of the Anderson model on 
the Bethe lattice was studied in detail in H. The criti- 
cal behavior found there is completely analogous to that 
obtained earlier for the Bethe lattice version of the a- 
models [Q,|l(J which can be derived from Wegner's n- 
orbital model with n 3> 1 states per lattice site O] 
(or, equivalently, considering a system of weakly cou- 
pled metallic granules 12 1). Exactly the same type of 
the transition and critical behavior on the Bethe lattice 
has been found (13] in the "toy" version of the supersym- 
metric er-model introduced in J14| , which is much simpler 
technically than the "true" cr-models of localization. In 
Ref. [pj| fluctuations of the local density of states (LDOS) 
near the Anderson transition were studied; as we will see, 
this question is closely related to the subject under dis- 
cussion. 

The aim of this paper is to study the model derived by 
AGKL via the supersymmetry technique. In fact, since 
many results for the critical behavior on the Bethe lat- 
tice are already known, we have just to apply them to 
the present model, properly identifying parameters and 
quantities of interest. There are, however, new features, 
which require an additional investigation. First of all, 
the model of AGKL has a large connectivity (branching 
number) m 3> 1. Secondly, we will be interested in the 
structure of the local spectral function in energy space; 
the question, which was not addressed in Refs. [p| ]l5[ . 

The outline of the paper is as follows. In SecH we 
formulate the model as derived by AGKL and identify 
the parameters and quantities to be studied. Sec. Ill is 



devoted to calculation of the average shape of the local 
spectral function on a Bethe lattice with large branch- 
ing number. In Sec.[TV| we evaluate the magnitude of the 
LDOS fluctuations in the region where these fluctuations 
are weak. In SecM we study the correlations of ampli- 
tudes of eigenstates on a tree-like lattice in the critical 
region. This information is used in Sec. VI, where we 



study fluctuations of the local spectral function and its 
typical shape in the vicinity of the Anderson transition, 
both in the phases of localized and extended states. In 



Sec. VII we discuss application of this results to the prob- 
lem of quasiparticle excitations in quantum dots. Our 



results are summarized in Sec. VIII. 



II. FORMULATION OF THE PROBLEM. 

For the sake of completeness and clarity, we recall the 
derivation of the model by AGKL ||. To describe the 
system of interacting fermions one should study eigen- 
states | a) of a generic many-body Hamiltonian of N pair- 
wise interacting particles : 
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where £ p are energies of the corresponding single particle 
states. When one neglects the interaction V the complete 
set of eigenstates is provided by Slater determinants ob- 
tained from the filled Fermi sea |0jv) as 
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As a result of interaction, the exact eigenstates \a) for 
interacting particles are superpositions of the Slater de- 
terminants | a) = £ M A£, } |tf# } ). 

To analyze the form of the line for single-particle exci- 
tations AGKL suggested to consider the hierarchy of the 
many-body states in the Fock space. Let us consider the 
system with N particles in its Hartree-Fock ground state 
|0jv) - The first level of hierarchy is formed by states 
with one particle: |l p ) = aJjOjv), second by the states 
with two particles and one hole (which we will denote as 
three-particle states for brevity): |3 9pr ) = aja p a r |0at), 
the third one by the states with three particles and two 
holes (to be denoted as five-particle states) and so forth. 
Neither of these states ( which we denote generally as \i) 
) is an eigenstate of Ti. Therefore if an extra particle is 
added to the system in such a way that a single-particle 
state |1) on the first level of hierarchy is formed, this state 
will decay to the states on the second level of hierarchy, 
which in turn may spread further in the Fock space. The 
process of the spread may be looked at as a quantum dif- 
fusion of a fictitious particle on a graph whose vertices i 
are formed by the states \i) at different levels of hierarchy. 
As such, it can be described by an effective one-particle 
tight-binding model characterized by the Hamiltonian: 
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with "site energies" e, = (i\H.\i) and hopping constants 

Vij = {i\n\j). 

Let us note that a problem of decay of an isolated level 
to a dense band in a situation when the level is directly 
coupled only to a small fraction of the states of the con- 
tinuum was studied by Akulin and Dykhne H . Though 
the model considered in |l6| is somewhat different from 



that studied in the present paper, many physical features 
are actually similar for both models. In particular, it was 
stressed in 16 that the dynamics of level decay can by 
affected by the phenomenon of Anderson localization. 

The spectral shape of a quasiparticle peak is obtained 
by projecting the initial single-particle state |1) onto ex- 
act eigenstates \a) ||: 



Pl (E)=^2\(a\l)\ 2 S(E-E 



(4) 



If we now consider the state |1) as a site of the tight- 
binding model, we see that p\ [E) is nothing else but the 
local density of states (LDOS) in the point 1 at energy 
E. i.e. the Fourier-transform 
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of the overlap of the spreading wave function with the 
initial state |1). In another context this quantity is called 
the strength function and is frequently used in nuclear 
and atomic physics to characterize spectra of complex 
systems, see e.g. |L7J for recent applications and further 
information. 

The interaction matrix elements couple the states of 
each generation to those of the preceding and of the 
following generations. In particular, a single-particle 
state is connected with three-particles states, but not 
with five-, seven-, . . .-particle ones. The density of one- 
particle states is equal to v\ — 1/A. Since we are in- 
terested in the regime, where one-particle states can be 
resolved or, in other words, the width is less than A, we 
can keep only the states with energies within a window 
(E — A/2, E + A/2) containing just one single-particle 
state. The density of three-particle states at given energy 
E S> A is given by 
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the factor 1/2 taking care of the identity of the quasi- 
particles. The ratio m = v$(E)/v\ = E 2 /4A 2 defines 
the effective branching number of the problem and is 
considered as a big parameter J18| . The important fact 
observed by AGKL is that this branching number stays 
parametrically the same for the following levels of hier- 
archy. For example, the density of the five-particle states 
vs(E) ~ E 4 /A 5 and each of these states is connected 
with of order of one three-particle state (in the chosen en- 
ergy window of the width A), whereas each three-particle 
state is connected to ~ E 2 /A 2 ~ m five-particle states. 
Following AGKL, we will neglect the change in numerical 
coefficients, and consider a tree structure with constant 
branching number m^> 1 (Bethe lattice) [|19[ . 

The Hamiltonian of the Anderson model on the Bethe 
lattice is determined essentially by the following three 
parameters: the branching number m, the characteristic 
magnitude V of the hopping matrix elements, and the 



width W of the distribution of the random site energies 
ej. They are related to the parameters of the original 
(quantum dot) problem as follows. The connectivity m 
is determined by the proliferation of the density of states 
from one level of the hierarchy to the next one, as dis- 
cussed above: 



m ~ v 3 (E)/ Vl ~ v 5 (E)/v 3 {E) 



E 2 /A 2 (7) 



The characteristic value of the hopping matrix element 
is equal to OH 



V ~ A 2 /E T ~ A/g 



(8) 



Finally, W can be determined by equating the mean spac- 
ing of the states in the next level of the hierarchy con- 
nected to a given state, l/v 3 {E) ~ A 3 /E 2 , to the corre- 
sponding quantity in the Anderson model on the Bethe 
lattice, W/m. This yields: 



W ~ mA 3 /E 2 ~ A 



(9) 



In order to check that restriction to the states in the 
energy interval A is justified, we can modify the model 
by including all states from a broader energy window 
W n = nA, with 1 <C n <C E/A. Then each three-particle 
state will be coupled to n one-particle states; each five- 
particle state to of order of n three-particle states etc. 
This modified model thus resembles the n-orbital model 
on the Bethe lattice. It has branching number given by 
the same expression (Q) and the coupling constant [ |ll|Jl2| 
a ~ (nV/W n ) 2 = (V/A) 2 ~ 1/g 2 , independent of n, 
which confirms that n is indeed an irrelevant parameter. 
Equivalence of the n-orbital and the Anderson (n = 1) 
models at a < 1 will be discussed in more detail else- 
where [ pl| . 

It is well known 0-pl that the Anderson model on the 
Bethe lattice exhibits the Anderson localization transi- 
tion when one of control parameters W, V, m varies. In 
our case, the change in the quasiparticle excitation energy 
E implies, see eq.(M) , the change of the branching num- 
ber m. Such a phase transition is reflected in a drastic 
change of statistical properties of eigenfunctions. In our 
preceding publications P,fl5[ we have discussed the dis- 
tribution function of LDOS and its typical spatial shape 
(for given energy) near the Anderson transition. Now, 
our interest is concentrated on a typical shape of p\ (E) 
in energy space, for a given "point" |1). 



III. AVERAGE LOCAL SPECTRAL FUNCTION. 

In this section, we study the average form of the local 
spectral density pi(E) in a site |1) of a Bethe lattice with 
large connectivity (branching number) m ^$> 1. We con- 
sider for definiteness a model with site-diagonal disorder 
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The averaging is performed over the random energies ej 
in all sites except the site |1), where the spectral den- 
sity is studied. Thus, we calculate the average form of 
the LDOS in a site with given value t\ of the random 
potential. The same object was studied previously in a 
random banded matrix model with strongly fluctuating 
diagonal elements pa). We will compare the results for 
the present case to those of the Ref. |28| in the end of the 
section. 

The average value of the LDOS is given by: M 



Pi 



(E) = -Re J dRRF m+1 {R 2 )e i{E ~ €l)R2/2 , (11) 



where F(R 2 ) is the solution to the following non-linear 
integral equation: 

F{S 2 ) = 1-S f dR f de 
x 1 {e)e l{E ~' )R2/2 F m {R 2 )VJ l {VRS) (12) 

Here 7(e) is the distribution of the random site energies 
ej, and J\(x) is the Bessel function. We look for an ap- 
proximate solution of eq.(12) in a Gaussian form 

F{R 2 ) = exp{- go i? 2 } (13) 

Substituting eq.jl^) in the r.h.s. of cq. (|l2|) , we get 

r.h.s.(12)= /de7(e)exp(-V 2 S* 2 j- — - 

= expL±V 2 S 2 Jde 1 (e) 
x[l + 9l S A + . . .} , 
where 



mgo — i(E — e)/2 
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mg -i(E-e)/2\ J 
From eqs.(12), (]13|) , and ( |l4| ) we get an equation for go 
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The Ansatz (O) is justified if the 0(S 4 ) correction in 
eq.dlj) can be neglected. Since 



F m {R 2 ) = exp{-m ff0 i? 2 }(l + m gi R 4 



(17) 



the integrals in the r.h.s. of eqs.dO), (O) are dominated 
by the region R 2 ~ 1/mgo- Furthermore, according to 
eq.dlS), gi ~ g 2 . Thus, the R 4 term in eq.(17h leads to a 
relative correction of order of 1/m, which is small in the 
case of large connectivity, m ~^> 1. Therefore, eq.(O) is 
indeed the proper approximation for m^$> I. 



Now we will analyze eq.(|l6|) in the two regimes: (i) 
(V/W) 2 < 1/m, and (ii) (V/W) 2 > 1/m, where W is a 
typical magnitude of the random energies e,. In the case 
(i) one can neglect go in denominator of §IQ) . This yields 



Re ff0 = ^(E) = - 
lmg = ±V 2 [dey(e)^- 



u 
2 



(18) 



Since a typical scale for 7(i?) is -f(E) ~ PF -1 we find 
.go ~ V^ 2 /!^. We see that under the condition (V/W) 2 <§; 
1/m, the quantity rng p can be indeed neglected in the 
denominator of eq.([L6|) as compared to£-t~ W. Then 
the LDOS, eq.([JT), takes the Breit-Wigner form 
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(we neglected the difference between m + 1 and to in the 
second line), i.e. it is a Lorentzian with a width 

T/2 = mi> = m7rV r2 7( J B) - mV 2 /W < W (20) 

The center of this Lorentzian is determined by the local 
random energy e\ (slightly shifted by u ~ Y). The width 
r is exactly the quasiparticle decay rate one would get 
by applying the Fermi Golden Rule: it is equal to 2ir 
times matrix element V squared times density of states 
into which the particle decays (m"/(ei) ~ m'y(E)). 

In the opposite case, (V/W) 2 ^> 1/m, we can neglect 
e in the denominator of eq. (Ffq) , which reduces then to 
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We find therefore 
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and 
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Thus, the LDOS in this limit has a semicircular form. 
The width of the semicircle is of order of JmV ^S> W jus- 



tifying omission of e < W in the denominator of eq.(16). 
The crossover from the Lorentzian to the semicircular 
shape of local DOS is completely analogous to that found 
for the model of random banded matrices with strongly 
fluctuating diagonal elements [E8[. The physical reason 
for this crossover is the same in both cases. When the 
disorder is sufficiently strong, an eigenstate is not spread 
uniformly over all the sites of the lattice, but rather is 



concentrated on a small fraction of lattice, formed by 
the sites which satisfy a kind of resonance condition with 
a given site |1). This leads to a Lorentzian form of the 
LDOS with a width T much less than the width W of the 
distribution of diagonal matrix elements (random ener- 
gies). Position of this Lorentzian for a site \i) is deter- 
mined essentially by the random energy a. The weaker 
is the disorder the larger is T. The crossover point is 
determined by the condition T ~ W. For a still weaker 
disorder the LDOS acquires the semicircular form, one 
and the same for all sites of the lattice. In this regime 
eigenstates are spread over all the lattice like in the Gaus- 
sian ensemble. 

Finally, we note that when the Bethe lattice model 
is used to describe excitations in a quantum dot only 
the first (Lorentzian) regime is relevant. Indeed, in this 
case one is interested in the regime T < A when the 
single-particle levels are resolved, whereas the energies 
€i of many-particle states included in the Bethe lattice 
model are distributed in the interval of the order of A, see 
Eq.(g). The Lorentzian form of the local spectral func- 
tion holds, of course, also for higher energies, E > Er, 
where the discreteness of the spectrum is unimportant 
and the Fermi Golden Rule is certainly applicable. To 
see this in our approach, one has to include in the model 
the states from a broader energy interval nA > T, as 
discussed in the end of SecO. However, we do not con- 
sider this region, where the single-particle states are not 
resolved and all the results (say, for the lifetimes) can be 
simply obtained in perturbation theory. 



IV. FLUCTUATIONS OF LDOS: REGIME OF 
WEAK FLUCTUATIONS 

In the preceding section, we have calculated the local 
spectral function pi (E) at fixed value of the random po- 
tential ei in a given site |1) but averaged over all the other 
random energies ej, i ^ 1. The main aim of the remain- 
ing part of the paper is to find out how a typical (rather 
than averaged) LDOS looks like. For this purpose, we 
study fluctuations of p\(E) in various regimes. It turns 
out that the parameter which distinguishes between the 
regimes (at fixed value of the branching number m | |29| ) 
is 



a - (V/W) 2 



(24) 



Our consideration in this section will be similar to that 
of the preceding one. This will allow us to describe fluc- 
tuations of the LDOS in the regime where they are weak, 
namely a ^$> 1/m 2 . 

To study the fluctuations of the LDOS, it is not suffi- 
cient to consider the integral equation ( |12| ) for the func- 
tion F(S 2 ) determining the average LDOS. Instead a 
more complicated equation on a function F(S 2 ,S 2 ) of 
two variables has to be considered INI: 



F(Sf,S%) = 4 [RidRi f R 2 dR 2 J t) (VR 1 S 1 )J (VR 2 S 2 ) 

x wSm Fm{RlRl) I del(e) 



x exp[i(E - e)(Rl - RJ)/2] 



(25) 



Let us look again for its approximate solution in the 
Gaussian form: 

F(Sl Si) = exp [-\s\{u - iv) + l -Sl(u + iv) ) CUW 



Substituting ( |26| ) in the r.h.s. of eq.(P5|), we get 

r v 2 s 2 i 



r.h.s. (25) = de~f(e)expl 
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2i m{u + iv) — (E — e) 



-m(u — iv) + E — e 

(27) 



Expanding this expression up to the terms of the order 
of S 2 and 5| one finds that (v + iu)/2 = go satisfies the 
same equation ( fig) obtained in sec. Ill where the average 
LDOS was considered. 

The function i^SfjSf) determines the fluctuations of 
the one-site Green's function 

G 1 (E) = J2\(°<M\ 2 (E-E a +iO)- 1 



(imaginary part of which is equal to —Trpi(E)) via the 
following equation [[15| 

F m +\Rl R\) ex P {z(£ ~ ci)(i^ - R 2 2 )/2} 

/oc /*oo 

du / dvV{u, v) exp{m(iZ? - R 2 2 ) - v(RJ + RJ)} 

(28) 

Here u and v are the real and imaginary part of the in- 
verse Green's function, 



[Gx^r^u- 



(29) 



and V(u, v) is the joint probability distribution of u and 
v. 

Substituting eq.((2|) into @), we find 



P(m, v) = 6(u - u )6(v - v ) , 



(30) 



where (neglecting the difference between m + 1 and m as 
in Sec.g) 



wo = -^{E - ei -mu) 



i\) 



mv 



(31) 
(32) 



We see that eqs.(|3l|), (32) together with the definition 
(p9|) reproduce precisely the average LDOS found in the 



preceding section. Therefore, the Gaussian approxima- 
tion (|26| ) corresponds to absence of any LDOS fluctua- 
tions, as is clearly seen from eq.(|30|). To calculate correc- 
tions to the Gaussian approximation we make the next it- 
eration and expand eq.(E7T) up to 0(S i ) terms. Assuming 
that we are in the Lorentzian regime, (V/W) 2 <C 1/m, 
we get after some algebra 



F(Sf,Sl) = exp 
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■55f(- 

7T "f(E 

4 mv 
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v A slsl + 



(33) 



Substituting this in the l.h.s. of cq.u28[) and expanding 
both sides of this equation in a power series in (R 2 — R 2 ) 
and (R 2 +R 2 ), we find the average values (v) — i>o, (u) = 
uq (as given by eqs.(pl|), (J32J)), and the variances 

V 2 



f-x f-y n 4 "/(E-mu) 
var(w) = var(w) = —V 



v 



(34) 



The parameters uq and vq determine the position and 
the width of the averaged (Lorentzian) LDOS whereas 
cq.pj) describes its fluctuations. The relative strength 
of the fluctuations is characterized by the ratio 



var(-D) 7r 3 j(E — inu) 



1 



W 



4mirj(E)V mV 



(35) 



Thus, the LDOS fluctuations are weak when a 3> 1/m 2 . 
In the opposite case the fluctuations are strong and the 
Gaussian Ansatz is not a good approximation any longer. 
Correspondingly, the expansion around it employed in 
this section breaks down. 



V. EIGENFUNCTIONS CORRELATIONS NEAR 
THE LOCALIZATION THRESHOLD 

Now we want to study a typical structure of LDOS in 
the vicinity of the Anderson transition point. For this 
purpose, we will need, however, information concerning 
the correlations in amplitudes of two different eigenfunc- 
tions with energies close to the mobility edge in the same 
spatial point. We will not use the fact that to ^ 1 here 
and so our calculation will be valid for arbitrary branch- 
ing number. Moreover, we believe that the results of this 
section are of even more general validity and are not re- 
stricted to tree-like lattices (see the discussion in the end 
of the section) . 

The quantity we want to calculate here is the corre- 
lation function a(r, E, to) representing the overlap of the 
eigenfunctions, 

<r(r,E, U ) = A 2 i? 2 -V)<E I^MI^MI 2 

x 5{E - E a )6{E + u - E fi )) , (36) 

where R 2 {u) is the level correlation function, 



R(Lu) = A 2 (J25(E-E a )S(E- 



m 



(37) 
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We are going to study it for energies E, E+uj close to the 
mobility edge (in the phase of delocalized states) so that 
the correlation length £ is large. We will assume however 
that the system size is larger than the correlation length, 
so that the system is in the critical regime, but not ex- 
actly in the critical point. We will compare a(r,E,uj) 
to the analogous correlation function containing a single 
eigenfunction, 



ir{r,E)=A(£\*a(r)\ A 6(E-E a )) 



which defines the inverse participation ratio, 



P(E) = J2*(r,E) 



(38) 



(39) 



In fact, the Bethe lattice is not a proper model to 
study quantities like n(r, E) since these quantities are 
well-defined only for finite systems. However, a finite 
Bethe lattice has a peculiar feature: most of its sites 
are located at the boundary. As the result, the proper- 
ties may depend crucially on boundary conditions even 
in the thermodynamic limit. There exists, however, an- 
other model - that of sparse random matrices (SRM) - 
which describes a lattice having locally a tree-like struc- 
ture similar to that of the Bethe lattice but possessing 
also large-scale loops ensuring that all sites are essen- 
tially equivalent. As a consequence, this model is free 
from the problem of boundary conditions. At the same 
time it is known that the localization transition in the 
SRM model is equivalent to that on the Bethe lattice 
pOf . The inverse participation ration ( p9| ) was studied in 
the delocalized phase of the SRM model near the local- 
ization threshold in Ref. pi] . It was shown there that its 
critical behavior is determined by an exponentially large 
scale C(E) emerging in the non-compact sector of the 
effective supersymmetric model, 



P(E) oc N- l C{E) 

C(E) oc exp{const|£: - E c \~ 1/2 } 



(40) 
(41) 



The scale C(E) has a physical meaning of the "correla- 
tion volume" |L5j and determines also the critical behav- 
ior of conductivity in the delocalized phase. We extend 
now the calculation of Ref. |21j] to find also the correlation 
function (Pq). 

To calculate the overlap function defined in Eq.(|3q) 
we follow [E2| and use the identity relating 7r(r, E) and 
a(r, E, u>) to advanced and retarded Green functions 
G R > A (r, E) = Ea=i I'MOI 2 (E±M- E a )~h 



2n 2 



A- l ir(r, E)6(oj) + A' 2 R 2 (u>)cx(r, E, u>) 



(42) 



= Re [(G R {r, E)G A {r, E + u) - G R {r, E)G R (r, E + lu))] 



Let us consider for definiteness the ensemble of real sym- 
metric SRM, corresponding to systems with unbroken 
time-reversal invariance. For any site index r = 1,...,7V 
we introduce one eight-component supervector $t = 

I $^, $> A ) consisting of two four-component supervectors 

$1 = (0<r,bi,(?!) ( T,b2,Cj,-0 ( T,/J, where indices a = R,A 
and 6, / are used to label advanced-retarded and boson- 
fcrmion subspaces, respectively. The ensemble-averaged 
products (G a G a ') for RSM model in the limit N > 1 
can be extracted from the paper |20f| and the Appendix 
D of the paper E3] and is given by: 



StrQA ; 



G a {r,E)G a (r,E + u>)}= [ 1 - -6 a , 



■o,b\4>o,b\<P<T> ,b\4><r> M) gT exp I 



< / DQ 

) gT = J d*(...)exp 



( iitpujN 



-E&L& + mg T ($) 



(43) 



The function g T ($) = g (&T*T$; $+£$) satisfies the 
integral equation: 



gT(^) = {[h F {^L^)-l]) 



ilT 



(44) 



where hp(t) = J dze~ ltz h(z) is the Fourier-transform of 
the distribution of nonzero elements of the SRM. The 
8x8 supermatrices T satisfy the condition T^LT = L 
where L = diag(l, 1, 1, 1, — 1, — 1, 1, 1) and belong to a 
graded coset space whose explicit parametrization can 
be found in |gj,g5[. The supermatrices Q are expressed 
in terms of T as Q = T~ 1 AT. At last, the matrix A = 
diag(l, 1, 1, 1, — 1, — 1, — 1, —1), and the density of states 
p is expressed in terms of the solution of the equation 
Eq.(||) as p(E) = -2g 0x /(nB 2 ), where B 2 = / dzh(z)z 2 
and g 0x = dg (x, y)/dxL >y = ; x = $ f $, y = $tL$. 

When deriving Eq. ( p3| ) , evaluation of a functional inte- 
gral by the saddle-point method has been employed, see 
details in |20|]23] ]. An accurate consideration shows that 
such a procedure is legitimate as long as: i) the matrix 
size N (playing in our model the role of the volume) is 
large enough (much larger than the coefficient C(E) de- 
termining the size dependence of IPR, see above); and 
ii) the energy difference u> is small enough (much smaller 
than C~ 1 (E)). Though C{E) is exponentially large near 
the transition point, it depends on the energy E only, 
so that when we keep E fixed and increase the system 
size N, the number of levels in the interval C~ 1 (E) gets 
arbitrarily large, since the level spacing scales as 1/AT. 

Expanding both sides of Eq.(|44|) over 'J, one can ex- 
bi4><rM4>cr' ,bi(f>a' m) rr in terms 0I tne matrix Q 



press 

as 



'St 



■,bl < P<J,bl ( Pcr.bl c Pcr,bl l 



'IT 



4! 

^>R,bl'i ) RM'i ) A,bl4>A,bl) 



290,xxQb x biQb t bi + ga.xyQbtbt + 90,yy 



— ~TT [90,xx (QbibiQbxbi + ^QbibiQbibi) + 
go,xy (Qb^QbrfJ + 90,yy] , 



(45) 



where g 0tXX = d 2 g /dx 2 \ x ^ =0 ; g , yy 



= d 2 g /dy 2 \ x ,y=o] 
go,xy = d 2 g /dxdy\ x , y= Q, and B 4 = / dzh(z]z 
allows us to represent right-hand side of Eq. 
following form: 



This 

in the 



2tt j 



A-\(r, E)S(u) + A- 2 R 2 (ui)a(r, E, w) 

4 



Ba 



goxxR&dQbibxQbibi + ^QbibiQbib! 



(46) 



where 



\ n RR n RR , n AA n AA ' 



>,,= / dQ(...)exp(^2^StrQA 



The integrals over Q— matrices are the standard ones [ |24| 
yielding: 



2iA 



/ i?.A n AR \ 

W6i6i^6ibi?Q - ^( w + z0 )' 

(Qb 1 b 1 Qb 1 b 1 /Q = (Qb 1 b 1 Qb 1 b 1 ) t 



(47) 



where i?2 ( W /A) is the level correlation function in the 
Gaussian Orthogonal Ensemble. Substituting this in 
Eq. (pi), we finally find 



a(r,E,oj) = -7r(r,J5) 



1 4g 0;a;x 
N 2 ir 2 p 2 B 4 



(48) 



The coefficient 1/3 in Eq.(f48D corresponds to the case 
of unbroken time reversal symmetry (orthogonal ensem- 
ble). For the unitary ensemble (broken time reversal sym- 
metry) the same consideration yields the coefficient 1/2 
instead, so that the general relation reads 



a{r,E,u) = 







+ 2 



n(r,E) , 



(49) 



where /3 is the conventional symmetry parameter equal 
to (3 — I (2) for the orthogonal (resp. unitary) ensem- 
bles. This relation between the overlap of two different 
eigenfunctions <j(r, E, to) and self-overlap n(r, E) is valid 
everywhere in the phase of extended eigenstates, up to 
the mobility edge E — E c , provided the number of sites 
(the system volume) exceeds the correlation volume. In 
particular, it is valid in the critical region \E — E c \ <C E c , 
where a typical eigenfunction is very sparse and self- 
overlap (hence, IPR) grows like exp (const|-E — i? c | -1 / 2 ) 

II- 

Eq.(B9|) implies the following structure of eigenfunc- 
tions within an energy interval SE = u) < C^ 1 (E). 



Each eigenstate can be represented as a product ^ (r) = 
ipi(r)<&E(r)' The function ^e{t) is an eigenfunction en- 
velope of "bumps and dips" which is smooth on a mi- 
croscopic scale comparable with lattice constant. It is 
the same for all eigenstates around energy E, reflects the 
underlying gross (multifractal) spatial structure and gov- 
erns the divergence of self-overlap at critical point. In 
contrast, ipi(r) is Gaussian white-noise component fluc- 
tuating in space on the scale of lattice constant. It fills in 
the "smooth" component $s(r) in an individual way for 
each eigenfunction, but is not critical, i.e. is not sensitive 
to the vicinity of the Anderson transition. These Gaus- 
sian fluctuations are responsible for the factor (3/(f3 + 2) 
(which is the same as in the corresponding Gaussian En- 
semble) in Eq.(|||). 

As was already mentioned, this picture is valid in 
the energy window SE ~ C~ 1 (E) around the energy 
E; the number of levels in this window being large as 
SE/A ~ NC^ 1 (E) 3> 1 in the thermodynamic limit 
N — ► oo. These states form a kind of Gaussian Ensem- 
ble on a spatially non-uniform (multifractal for E — > E c ) 
background $g(r). Since the eig enfunction correlations 
are described by the formula (Hfl) , which has exactly the 
same form as in the Gaussian Ensemble, it is not surpris- 
ing that the level statistics has the WD form everywhere 
in the extended phase [M . 

Let us make a side remark here. We believe on physi- 
cal grounds that the same picture should hold for a con- 
ventional d-dimcnsional conductor. First of all, the gen- 
eral mechanism of the transition is the same in d < oo 
and d — oo models. Furthermore, the sparsity (multi- 
fractality) of eigenstates near the transition point takes 
its extreme form for d = oo models |Lq |, so that since 
the strong correlations (|49]) take place at d = oo it 
would be very surprising if they do not hold at finite 
d as well. Finally, Eq.([49|) was proven by an explicit 
calculation in the weak localization regime |22|, where 
a(r,E,uj) = ^n{r,E) = V"~ 2 [l + U(r, r)], with V be- 
ing the system volume and II(r, r) the diffusion propaga- 
tor. Therefore, we believe that also for a (i-dimensional 
conductor the WD statistics applies everywhere in the 
delocalized phase, if the system is large enough. Fur- 
thermore, this implies that exactly in the critical point 
the level repulsion has a conventional uj^ form on a scale 
u> ~ A. We refer the reader to a separate publication 
p7J , where these questions are discussed in more detail. 



VI. FLUCTUATIONS OF THE LDOS: CRITICAL 
REGION. 

Now we are prepared to study the structure of the 
LDOS in the vicinity of the localization transition. We 
note first that the parameter a defined by eq. (E4) is noth- 
ing else but the coupling constant of the corresponding 
tr-model. In fact, to derive a-model rigorously one has to 
introduce n 3> 1 states per lattice site |ll],|l^], whereas 



the Anderson model corresponds to n = 1. However, the 
both models have essentially equivalent critical behav- 
ior, see below. The cr-model is defined by the action (we 
consider the unitary symmetry case for definiteness): 



s{Q} = ^J2 StT Q*Qj + £ J2 StrA & ' 



(50) 



(ij) 



where the summation in the first term goes over the pairs 
of nearest neighbor lattice sites, Qi — T^~ ATj is a 4 x 4 
supermatrix belonging to certain coset space and satisfy- 
ing the constraint Q 2 = 1, e = — ittvuj /2, v is the density 
of states, A = diagjl, 1,-1, —1}, and Str denotes the su- 
pertrace. The reader is referred to Refs. |24|-|26|| for more 
detailed information on the supersymmetry formalism. 

While considering the vicinity of the Anderson tran- 
sition point we prefer to work within the a— model for- 
malism. The cr-model formulation has a considerable ad- 
vantage: its behavior is determined by a single coupling 
constant a. This should be contrasted with the original 
Anderson model controlled by energy E and the whole 
distribution function 7(e)). However, the mechanism un- 
derlying the transition and all the essential features of 
the critical behavior are exactly the same for both the 
cr-model and the Anderson model on the Bethe lattice, 
see Refs. f@-RM3 . 



When formulated on a Bethe lattice the nonlinear cr- 
model can be reduced to a certain non-linear integral 
equation for a function Y(Q) on the coset space defined 
as 



Y(Qi) = Jl[DQ i 



-S{Q} 



(51) 



Because of the invariance properties the function Y(Q) 
depends only on the two scalar variables (eigenvalues) 
1 < Xi < 00 and — 1 < A 2 < 1. Moreover, in the lo- 
calized phase (as well as close to the transition point in 
the delocalizcd phase) only the dependence on the "non- 
compact" variable Ai is essential. In particular, in the 
localized phase the integral equation Eq. (pTj) in the limit 
rj — > takes the form 

/oo 
dt'L a {t - t') exp{-2e*'} 2 / m (t') , (52) 
-OO 

where t = s\\ and 

1 lo 

L a (t) = (JL) e */2 e -acosht [ sinhacoshi 



I cosh a 



sinha 
2a 



(53) 



In the localized regime (a < a c ), eq.(|52|) has a solution 
in the form of a kink with the asymptotics y(t) ~ 1 at 
t — » —00 and y(t) ~ at t — ► +00. In contrast, in the 
delocalized phase (a < a c ), eq.(|52]) has a trivial solution 
y(t) = only. The Anderson transition point a c is thus 



determined by a condition of the stability of the kink 
solution, which has the form p,|lO|| 



mw ae (l/2) = l ; 
Wa ($) = f dte~ et L a (t) 



(54) 
(55) 



For the nonlinear cr-model with orthogonal symmetry 
(corresponding to the case of unbroken time reversal in- 
variance) the underlying structure is completely analo- 
gous, except for a more complicated form of the function 
L a (t) P,p0[. However, all the gross properties and the 
results are essentially the same as in the unitary case. 

The technically simplest (but still non-trivial) "toy" 
variant of the supersymmetric nonlinear cr-model was 
introduced by Zirnbauer in |14| and called "hyperbolic 
superplane" (HSP). When formulated on the Bethe lat- 
tice it has the critical behavior completely equivalent to 
that of the true nonlinear cr-models of Anderson local- 
ization p3[ . In particular, the "invariant" phase (which 
should be associated with the localized states regime) is 
described by the same integral equation (|2|), with the 
kernel L a (t) given by 



(56) 



La ^ = ^n? — rr exp ( o ~ acoshi 

2K 1 / 2 (a) \2 



The kernel (|36j) has the same gross features as (j53 
Moreover, we are able to show that the position of the 
mobility edge in the Anderson model on the large connec- 
tivity Bethe lattice is exactly given by Eqs.(^4|), (55) with 
the same kernel Eq.(56). Derivation of this result, and 
more generally, demonstration of the equivalence of the 
a— model and the Anderson model in the regime a«l 
will be published elsewhere f3l]| . Let us stress, that this 
is exactly the re gion , where the transition for m ^$> 1 
happens, see Eq.( p8|) b elow. 

Substituting eq.(|56j) into eqs.(|54]), (|55|), one gets the 
following equation for the transition point [13]: 



, K (a c ) 

K 1/2 (a c ) 



1 



(57) 



where K v {z) is the modified Bessel function. Using now 
■^\/%{ z ) = (7r/2z) 1 / 2 e _z and the small-z behavior of 
K (z), K (z) ~ -log(z/2), we find 



8m 2 log m 



to > 1 



(58) 



Let us first analyze the structure of the LDOS pi(E) 
in the phase of localized states, a < a c . In this regime, 
the spectrum is discrete, so that pi(E) is given by a sum 
of (5-functions with positions in E = E a and residues 
|(a|l)| 2 . Only those eigenstates |a), which are localized 
close to the site |1), give appreciable contributions to this 
sum. The effective number of such states is measured by 
P _1 (ei), where P(E) is the inverse participation ratio, 



P(i5) = X>|1> 



"£\{a\l)\*6(E a -E) 



so that 



(59) 



Rewriting P(E) in the second form (more convenient for 
an analytical study), we used the fact that the states \a) 
contributing appreciably to p± (E) have energies E a close 
to e\. The latter statement follows immediately from the 
results of Sec. [II, where the average LDOS was shown to 
be a narrow Lorentzian around ei, see eq. (|l9|) . 

The IPR takes values between 1 and 1/N correspond- 
ing to complete localization and delocalization of eigen- 
functions respectively. In the localized phase, P(E) is 
finite, whereas in the delocalized one it is proportional to 
1/N (inverse system volume), and thus is equal to zero 
in the thermodynamic limit. In a d- dimensional conduc- 
tor, when E approaches the transition point from the 
phase of localized states, P(E) vanishes continuously: 
P{E) oc \E — E c \* 2 , where tt2 is one of the set of the ex- 
ponents determining the multifractal structure of eigen- 
functions |3^]. In contrast, on a Bethe lattice this hap- 
pens in an abrupt manner [p|- |l0| , [30"| ] : P(E) has a finite 
limiting value P c as E approaches E c from the side of 
localized states (a — > a c — in terms of the er-modcl) 
and then jumps discontinuously to zero in the extended 
phase (a — a c + 0). Therefore, the effective number of 
5-peaks in p\ (E) in the localized phase stays finite as we 
approach the Anderson transition point. 

This result was obtained in Refs. [ra-flOpQ], where the 
lattice connectivity was considered as a number of order 
of unity. Let us study now what happens with the lim- 
iting value P c for large connectivity m ^> 1. The IPR is 
given by the following expression JTo|,|30[ 



P 



dte t exp{-2e t } 2 / m+i (t) 



(60) 



where y(t) is the solution of eq. (|5^) . For m ^> 1 the 
transition happens at a c <C 1, see eq.(p8|). It is easy to 
see that the function L a (t), eq. (p6[) , has in this regime a 
sharp maximum around t = log(l/a) 3> 1 with a width of 
the order of unity. Taking into account the normalization 
property, J dtL a (t) = 1, we have thus L a (t) « 6(t — 
log(l/a)). Substituting this in eq. (|52|) , we find that the 
factor exp(— 2e* ) in the r.h.s. implies that the kink in 
the function y(t) is located at t ~ log(l/a). We know 
from the analysis of Refs. g [bj,[]j|,|30| that the small-* 
behavior of y(t) near the transition point is the following: 



Vc(t) 



1 



,(t-t )/2 



t < t , 



(61) 



where to is a constant determining the position of the 
kink. Wc find thus e*° ~ l/a c , up to a factor of the 
order of unity. Substituting this in cq.(BG), we get 



P 



d£e*exp{-2e'} 



1 



,(t-*o)/2 



(62) 



1 - P c = 2me- to/2 / dtexp{3t/2 - 2e 1 } 



-to/2 



.1/2 



logTTJ 



(63) 



Therefore, at m ^$> 1 not only IPR has a non-zero limit 
at the transition point but this limiting value is close to 
unity, i.e. to the upper bound corresponding to extremely 
localized states concentrated on a single site. Roughly 
speaking, the system undergoes a transition directly from 
the deeply localized phase (P close to unity) to the ex- 
tended phase (P = 0). 

The above results, eqs. (pl|)— (|63|) can be also repro- 
duced by solving iteratively the self-consistency equation 
(p2]). Numerical results obtained in Refs. |l(j,|3(| show a 
clear tendency of increase of P c with m, in full agreement 
with eq.(|63|). Therefore, when the system approaches 
the Anderson transition point from the localized side, 
the strength function pi (E) has a form of the sum of 5- 
like peaks with almost the whole of the spectral weight 
(except a small fraction of the order of 1/logm) concen- 
trated in one peak. 

When the system is driven through the critical point 
into the delocalized phase, this picture evolves gradually. 
Namely, the peaks get broadened with a width depending 
on the distance to the critical point and vanishing (in the 
thermodynamic limit) at the critical point. It is clear 
that the width is given by the scale C~ 1 (E), within which 
the eigenfunctions are fully correlated (see secM). This 
scale is exponentially small in the vicinity of the critical 



point, see eq 
constant a [|[ 



C 1 (a) oc exp • 




or in terms of the a-model coupling 



-ci 



Of 



-1/2- 



a > a r 



(64) 



Such a critical behavior is a peculiar feature of the tree- 
like lattices, and is replaced by a conventional power- 
law critical behavior for a d-dimensional systems with 
d < oo [|15|]. Therefore, in the critical vicinity of the 
transition point (on the delocalized side), the peaks in 
p\ (E) have an exponentially small width of the order of 

(a), eq.(p4). This is also confirmed by the behavior 

j) of the IPR in this region. 

Let us study now the behavior of the coefficient c\ of 
the critical behavior (J6J) in the limit of large connectivity 
m ^> 1. This can be again the most easily done for the 
HSP model. The condition on the transition point is 
given by eq. (|5^) , and the scale C(a) is given by (l3| 



C(a) 



o*h 



where r\ is determined from the equation 



m 



Kjr, (a) 
#1/2 (a) 



(65) 



(66) 



Expanding eq. fl6q ) in i] and a — a c and using the asymp- 
totic behavior of the modified Bessel function, 



K iri (a) ~ - 



1(^(2) + r H^)(2 



a< 1 



r(«j) = — r(i + i?7)~— (1-iTjC) , 77 <i 

(C being the Euler's constant), we reduce eq.(|66|) to the 
following form: 



(a - a c ) (67) 




Using now the formula (pq) for a c , we find 



6 , 

— m[a — a c ) 



1/2 



(68) 



so that according to eq.(|65|), 




(69) 



Therefore, the coefficient ci in eq.(p4|) behaves as c\ ~ 
C2logm, with a factor ci of the order of unity. Ana- 
lyzing the derivation of eq. (BS ) , we find that the critical 
behavior in the form ( |69| ) , (p4 ) is valid for a — a c <C a c , 
where the fluctuations are exponentially strong. At 
a — a c ~ a c the quantity C(a) as given by eq. (p4|) is 
still large, C(2a c ) ~ e Cl = m C2 , in view of m ~^> 1. How- 
ever, C(a!) is not exponentially large anymore and ceases 
to be the leading factor determining the critical behavior 
of most of the relevant quantities, such as the LDOS mo- 
ments and the conductivity. For this reason, it turns out 
to be difficult to extend the above consideration onto the 
> a c . We know, however, from Sec. IV that 



region a — a c 

the region of relatively strong LDOS fluctuations extends 

up to a ~ 1/ra 2 ~ a c \ogm. 



VII. APPLICATION TO THE PROBLEM OF A 

QUASIPARTICLE LINE SHAPE IN QUANTUM 

DOTS 

Let us now translate the obtainter results into the con- 
text of the problem of the one-particle excitations in a 
quantum dot. This can be straightforwardly done by 
using the relations (J7]), (g), @, and (gj) between the 
parameters of the two problems. The localization transi- 
tion then corresponds to the energy E c ~ A(g/\ogg) 1 ^ 2 
H. In the localized region (excitation energy E below 
E c ) the bare states forming the basis of the Fock space 
mix only weakly to each otner, so that the exact eigen- 
states are close to the bare ones. In particular, admixture 



of many-particles states to a single particle one is weak. 
Therefore, only one exact eigenstate will contribute es- 
sentially to the spectral decomposition (Eh of a single- 
particle state, see fig la. In the delocalized domain (E 
above E c ) exact eigenstates are superpositions of many 
bare ones. This is in full analogy with a delocalized state 
of the tight-binding model, which covers many (infinitely 
many in the thermodynamic limit) sites of the lattice. As 
a result, there are many exact eigenstates contributing to 
the strength function (|J) . The corresponding envelope is 
of irregular (strongly fluctuating) shape in the intermedi- 
ate (critical) regime E c < E < E' c , with E' c ~ Ag 1 / 2 , and 
acquires a Breit-Wigner form at E ^> E' c . The width of 
this Breit-Wigner envelope (spreading width of the one- 
particle state) is given by eq.(|20|) (Golden Rule), yielding 
r~ A{E/Ag) 2 . 

Finally, let us remind that what we considered 
throughout the paper was a tree-like model with a con- 
stant coordination number. In reality, however, the num- 
ber of states to which a state of the n-th generation is cou- 
pled, decreases with increase of n. Let us briefly discuss 
how this is expected to modify the results. The transi- 
tion will be now smeared into a crossover, since the tran- 
sition point gets "generation-dependent" (see also recent 
preprint |33[). Furthermore, in the constant coordination 
number approximation all states on the energy shell (i.e. 
with energies within the spreading width ~ V around E) 
get mixed at E 3> E' c . In contrast, now only first few gen- 
eration will get mixed under this condition (that will be 
however sufficient to produce the Breit-Wigner envelope 
of the spectral function (Bj)). Admixture of higher gen- 
erations will require higher energies. To estimate, when 
the complete mixing of the states on the energy shell hap- 
pens, we note that the density of states of the generation 
n is equal to 



v 2n +i{E) 



1 



E 



n\{n + l)\{2n)\A \A 



which is a direct generalization of eq.(J6|). Maximizing 
this expression, we find that a typical many-particle state 
belongs to a generation with n ~ (E/2A) 1 / 2 , with typ- 
ical energies of quasiparticles ~ E/n ~ {EA) 1 / 2 . The 
level spacing of the states to which this one is coupled 
is ~ A(A/_E) 3 / 2 . Comparing this to the typical value 
of the matrix element, V ~ A/g, we conclude that the 
full mixing of the states on the energy shell (ergodicity) 
will be reached at E > £"> witl1 E " ~ A3 2 / 3 . The same 
estimate for the "chaotization border" was obtained very 
recently by Jacquod and Shepelyansky p3. 



VIII. SUMMARY 

In this paper we have studied in detail the structure of 
the average and typical local density of states in a tight- 
binding model on a tree-like lattice with a large branch- 
ing number m ^> 1. In the framework of the mapping 
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recently suggested in Ref. || this local density of states 
describes the shape of the quasiparticle excitation line 
in a finite Fermi system (e.g., quantum dot). We have 
exploited the supersymmetry approach to the problem 
developed previously, see Refs. p-jlQ,OB0|, and some of 
the results obtained in these papers. 

The results depend on the relation between two di- 
mensionless parameters: the branching number to and 
the coupling constant a ~ (V/W) 2 (V being the typ- 
ical magnitude of the hopping matrix element in Fock 
space and W the width of the distribution of random 
site energies ej). The relation of these parameters to 
those of the original quantum dot model can be found 
in Sec.||. When a ^> 1/m the LDOS in a given site 
i averaged over the distribution of all random energies 
of other sites ej, j =/= i has a semicircular form. In 
the opposite case (relevant to the quantum dot model), 
a <C 1/ffi, the averaged LDOS has a Lorentzian form 
with a width given by the Golden Rule, see eq. (p0|) . How- 
ever, the typical LDOS is close to its average value only 
for a > 1/m 2 . In the opposite case, a <C 1/m 2 , the 
LDOS fluctuations are strong, and the averaged LDOS 
is not representative, see Fig.l. This is related to the 
existence of the Anderson localization transition at the 
point a = a c ~ l/(m 2 log m). 

On the insulating side of the transition, a < a c , the 
LDOS is given by a discrete sum of ^-function peaks. 
The effective number of such peaks is characterized by 
1/P, P being the inverse participation ratio. It is known 
H |bj,g0) that on the Bethe lattice the inverse participa- 
tion ratio has a finite limiting value P c when the system 
approaches the transition point from the localized phase. 
We have shown that in the limit to 3> 1 this limiting 
value is close to unity, 1 — P c ~ 1/ log to <C 1, so that 
almost all the spectral weight (0) (except a small part of 
the order of 1/ log to) is concentrated in a single 6-peak. 
Roughly speaking, the system undergoes a transition di- 
rectly from the deeply localized phase to the extended 
one. 

When the system is driven through the critical point 
into the phase of extended states the peaks get broad- 
ened, with their width depending on the distance to the 
critical point and vanishing (in the thermodynamic limit) 
at the critical point. The width is determined by the 
scale C^ 1 (E) such that for energy separations smaller 
than C~ l (E) different eigenfunctions are fully correlated 
(see sec.[V]). As a result, the width is exponentially small 
in the vicinity of the critical point, see eqs.(fll|), (64). 
This is also confirmed by the behavior (MO) of the inverse 
participation ratio in this region. 

Our results by and large confirm the picture presented 
by AGKL ||. We have, however, quantified many fea- 
tures of the problem by using the supersymmetry ap- 
proach and some of the results obtained earlier in the 
framework of this method. 

We disagree with AGKL concerning their statement 
that in the delocalized phase the level correlation func- 
tion of the tree-like model with fixed coordination num- 



ber is not of the Wigner-Dyson form close to the tran- 
sition point because exact eigenstates are very sparse 
and "do not talk to each other" fa]. We have shown 
in Sec.|V| that the eigenfunctions of the sparse random 
matrix model which are close in energy are strongly cor- 
related in this regime. This explains the Wigner-Dyson 
form of the level correlation function proven in one of our 
earlier publications |20[| . 

For the quantum dot problem the obtained results im- 
ply the Fock-space derealization of a single-particle ex- 
citation at excitation energies E > E c ~ A^g/logg) 1 / 2 
and formation of regular-shaped Breit-Wigner envelopes 
at E > E' c ~ Ag 1 / 2 . Taking into account decrease of 
the coordination number in higher generations allows to 
estimate an energy, at which the complete ergodicity on 
the energy shell is restored as E" ~ Ag 2 / 3 . 

Finally, we would like to mention that an extensive 
study of statistical properties of eigenstates of a finite 
system of interacting fermions was undertaken recently 
by Flambaum, Izrailev, Casati, and Gribakin p5| , p6[ . 
In particular, these authors discuss crossover from the 
regime of strongly fluctuating local spectral density to 
that characterized by small fluctuations, which is similar 
to the questions addressed in the present paper fl37|. 
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